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ABSTRACT 

We  report  measurements  of  the  [C  n]  157.74  pm  fine-structure  line  in  a  sample  of  seven  ultraluminous  infrared 
galaxies  (ULIGs)  (LIR  >  10 12  L0)  with  the  Long  Wavelength  Spectrometer  on  the  Infrared  Space  Observatory. 
The  [C  ii]  line  is  an  important  coolant  in  galaxies  and  arises  in  interstellar  gas  exposed  to  far-ultraviolet  photons 
(hv  >  11.26  eV);  in  ULIGs,  this  radiation  stems  from  the  bursts  of  star  formation  and/or  from  the  active  galactic 
nuclei  that  power  the  tremendous  infrared  luminosity.  The  [C  n]  158  pm  line  is  detected  in  four  of  the  seven 
ULIGs;  the  absolute  line  flux  (about  a  few  times  1 0  20  W  cirT2)  represents  some  of  the  faintest  extragalactic 
[C  n]  emission  yet  observed.  Relative  to  the  far-infrared  continuum,  the  [C  ii]  flux  from  the  observed  ULIGs  is 
~10%  of  that  seen  from  nearby  normal  and  starburst  galaxies.  We  discuss  possible  causes  for  the  [C  n]  deficit, 
namely  (1)  self-absorbed  or  optically  thick  [C  ii]  emission,  (2)  saturation  of  the  [C  n]  emission  in  photodissociated 
gas  with  high  gas  density  n  (»  3  x  103  cirT3)  or  with  a  high  ratio  of  incident  UV  flux  G0  to  n  ( GJn  a  10  cm3), 
or  (3)  the  presence  of  a  soft  ultraviolet  radiation  field  caused,  for  example,  by  a  stellar  population  deficient  in 
massive  main-sequence  stars.  As  nearby  examples  of  colliding  galaxies,  ULIGs  may  resemble  high-redshift 
protogalaxies  in  both  morphology  and  spectral  behavior.  If  true,  the  suggested  [C  ii]  deficit  in  ULIGs  poses 
limitations  on  the  detection  rate  of  high-z  sources  and  on  the  usefulness  of  |C  n]  as  an  eventual  tracer  of 
protogalaxies. 

Subject  headings:  galaxies:  active  —  galaxies;  ISM  —  galaxies:  starburst  —  infrared:  galaxies  — 
infrared:  ISM:  lines  and  bands  —  ISM:  atoms 


1.  INTRODUCTION 

The  157.74  pm  2PV2  — ►  2Pl/2  fine-structure  line  of  C+  is  the 
single  brightest  emission  line  in  the  spectrum  of  most  galaxies, 
providing  as  much  as  1%  of  the  total  far-infrared  (FIR)  lu¬ 
minosity  (see,  e.g.,  Stacey  et  al.  1991b  and  references  therein). 
In  the  interstellar  medium  (ISM)  of  galaxies,  the  [C  ii]  158  pm 
line  traces  gas  exposed  to  stellar  far-ultraviolet  (far-UV)  pho¬ 
tons  with  energy  greater  than  1 1 .26  eV,  the  ionization  potential 
of  neutral  carbon.  In  such  regions,  atomic  hydrogen  and  elec¬ 
trons  collisionally  excite  ground-state  C+  ions,  producing 
[C  ii]  radiation  that  cools  the  gas.  The  [C  ii]  transition  has  a 
critical  density  ncr  of  3  x  103  and  50  cirT3  for  collisions  with 
hydrogen  and  electrons,  respectively  (Flower  &  Launay  1977; 
Hayes  &  Nussbaumer  1984),  and  modest  energy  requirements 
for  excitation,  i.e.,  7'levcl  ~  92  K  above  the  ground  state.  These 
attributes,  combined  with  a  high  abundance  of  interstellar  car¬ 
bon  (e.g.,  C/H  =  1.4  x  10~4;  Sofia  et  al.  1997),  make  the 
[C  ii]  line  the  dominant  coolant  in  UV-exposed  atomic  gas  over 
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a  wide  range  of  physical  conditions  (Hollenbach,  Takahashi, 
&  Tielens  1991;  Wolfire  et  al.  1995).  Such  environments  in¬ 
clude  much  of  the  warm  (IS  100  K),  dense  (n  >100  cnU3) 
atomic  ISM  between  H  ii  regions  and  molecular  clouds,  i.e., 
photodissociation  regions  (PDRs),  as  well  as  the  cold 
(T  ;£  50  K),  diffuse  phases  of  the  ISM.  The  [C  n]  line  can  also 
be  produced  in  the  warm  ( T  ~  8000  K)  neutral  ( njn  ~  0.1)  and 
ionized  (njn  ~  1)  diffuse  phases. 

As  examples  of  gas-rich  galaxies  characterized  by  strong  star 
formation  activity,  ultraluminous  infrared  galaxies  (ULIGs) 
(Lm  >  1012  L@)  seem  to  offer  the  prime  ingredients  for  bright 
[C  ii]  emission.  The  tremendous  infrared  (IR)  luminosity  of 
ULIGs  is  ~103  times  that  of  the  Milky  Way  Galaxy  and  is 
thought  to  be  due  to  enhanced  star  formation  (“starburst”) 
and/or  an  active  galactic  nucleus  (AGN),  which  produces  UV 
radiation  that  is  absorbed  by  dust  and  reradiated  in  the  IR.  The 
starburst  and/or  AGN  form  as  a  result  of  the  merger  of  two 
gas-rich  spiral  galaxies,  thereby  creating  a  large  concentration 
of  molecular  gas  (MH,  >  101()  M0)  in  the  core  which  fuels  the 
starburst  and  AGN  (e.g.,  Norman  &  Scoville  1988).  The  star- 
burst/AGN  activity  likely  dominates  the  excitation  and  global 
energetics  of  all  phases  of  the  ISM  in  these  galaxies. 

To  compare  the  ISM  in  ULIGs  with  that  in  less  luminous 
galaxies,  we  have  acquired  spectra  of  the  [C  ii]  158  pm  line 
in  ULIGs  using  the  Long  Wavelength  Spectrometer  (LWS) 
(Clegg  et  al.  1996)  on  the  Infrared  Space  Observatory  (ISO) 
(Kessler  et  al.  1996).  In  addition  to  tracing  the  multiphase  ISM, 
the  FIR  [C  ii]  line,  compared  to  optical  lines,  for  example,  is 
much  less  affected  by  the  high  extinction  (Av  S  1000;  Sanders 
&  Mirabel  1996)  typical  of  the  heavily  obscured  nuclear  regions 
in  ULIGs.  The  observations  presented  here  constitute  part  of 
the  LWS  consortium  guaranteed-time  “Central  Program”  on 
IR-bright  galaxies.  In  a  future  paper,  we  will  discuss  the  results 
of  an  ongoing  ISO  study  that  examines  the  [C  n]  emission  in 
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TABLE  1 
Observing  Log 


Source  a2000  S2000  AOTa  Observation  Date 


Arp  299  .  11  28  32.3  +58  33  43.4  LWS01  1996  May  15 

IRAS  12071-0444  .  12  09  45.3  -05  01  13.3  LWS02  1996  Jul  19 

Mrk  231  .  12  56  14.2  +56  52  24.9  LWS01  1996  May  15 

IRAS  15206+3342  .  15  22  38.1  +33  31  33.3  LWS02  1997  Sep  5 

IRAS  15250+3609  .  15  26  59.3  +35  58  37.2  LWS02  1996  Sep  20 

Arp  220  .  15  34  57.2  +23  30  11.2  LWS01  1996  Aug  20 

NGC  6240  .  16  52  58.8  +02  24  04.3  LWS01  1996  Aug  21 

IRAS  22491-1808  .  22  51  49.3  -17  52  24.1  LWS02  1996  May  21 

IRAS  23365+3604  .  23  39  01.2  +36  21  10.0  LWS02  1997  Dec  7 


Note.— Units  of  right  ascension  are  hours,  minutes,  and  seconds,  and  units  of  declination 
are  degrees,  arcminutes,  and  arcseconds. 

*  AOT  =  Astronomical  Observation  Template.  LWS01  =  medium-resolution  (R  ~ 
200)  spectrum  covering  the  full  43-197  /xm  range  of  the  LWS;  LWS02  =  medium- 
resolution  line  spectrum,  in  this  case  with  a  3.4  /xm  width  and  centered  on  the  [C  ii]  line 
wavelength. 


an  added  sample  of  ULIGs  and  follows  up  on  our  initial  guar- 
anteed-time  program.  Here,  we  also  include  for  purposes  of 
comparison  ISO  [C  ii]  measurements  for  two  luminous  IR  gal¬ 
axies,  Arp  299  and  NGC  6240  (XIR  =  8  and  7  x  10"  L0; 
Sanders,  Scoville,  &  Soifer  1991),  which  are  examples  of  col¬ 
liding  galaxies  but  are  slightly  less  luminous  than  the  ULIGs. 

2.  OBSERVATIONS  AND  RESULTS 

In  Table  1,  we  summarize  the  observing  log  for  the  ULIGs, 
along  with  Arp  299  and  NGC  6240,  included  in  the  ISO  LWS 
consortium  Central  Program  on  IR-bright  galaxies.  Observa¬ 
tions  of  the  [C  ii]  158  /tm  fine-structure  line  in  each  of  the 
galaxies  used  the  ISO  LWS  medium-resolution  (AX  =  0.6  /xm 
at  158  /xm)  grating  mode  (Clegg  et  al.  1996).  Lor  the  brighter 
sources  (F100fim  >  25  Jy),  we  obtained  full  43-197  /xm  spectra, 
from  which  we  extracted  [C  ii]  spectra.  Lor  the  fainter  objects, 
we  acquired  spectra  in  the  vicinity  of  the  [C  ii]  line  only. 

We  co-added  processed  spectral  scans,  extracted  line  profiles, 
and  measured  line  fluxes  and  uncertainties  using  the  ISO  Spec¬ 
tral  Analysis  Package  (ISAP)  and  IRAL  Pipeline  processing 
of  LWS  spectra  performs  wavelength  and  flux  calibration,  re¬ 
moves  cosmic -ray  glitches,  subtracts  dark  current,  and  corrects 
for  instrumental  responsivity  variations  (Swinyard  et  al.  1996). 
In  Table  2,  we  list  the  measured  [C  ii]  line  fluxes  with  the 
statistical  errors  and  3  a  flux  upper  limits,  where,  in  all  cases, 
the  line  was  unresolved.  The  upper  limits  were  calculated  as¬ 
suming  a  Gaussian  line  with  the  effective  instrumental  profile 
and  a  3  a  amplitude.  When  detected,  the  [C  n]  line  appeared 
at  the  expected  wavelength  adjusted  for  redshift.  At  the  red- 
shifted  wavelengths  of  the  [C  n]  line,  the  [C  n]  measurements 


are  uncontaminated  by  Galactic  [C  ii]  emission  for  all  sources. 
Uncertainties  in  the  flux  measurements  were  dominated  by  flux 
calibration  errors  arising  from  uncorrected  fluctuations  of  de¬ 
tector  responsivities  with  time  and  inaccurate  dark  current  sub¬ 
traction.  These  errors  are  estimated  at  ±25%  based  on  over¬ 
lapping  data  in  neighboring  detectors  from  the  full  grating 
scans.  The  LWS  beam  size  is  69"  at  158  /xm  (Burgdorf  et  al. 
1997).  All  sources  in  this  study  were  spatially  unresolved,  with 
the  beam  encompassing  all  of  the  IR  emission.  Lor  each  source, 
we  list  in  Table  2  the  LIR  flux  FmR,  defined  here  as  the  emission 
from  42.5  to  122.5  /xm  as  measured  by  IRAS  and  computed 
according  to  the  formula  FF1R  =  1.26  x  1CT18  [2.58  Fu( 60  /xm) 
+  F/lOO  /xm)]  W  cm~2,  where  T„(60  /xm)  and  T„(100  /xm)  are 
the  flux  densities  in  janskys  at  60  and  100  /xm,  respectively 
(Helou  et  al.  1988). 

As  shown  in  Table  2,  the  [C  ii]  line  has  been  detected  in 
four  of  the  seven  ULIGs  observed  to  date,  plus  Arp  299  and 
NGC  6240.  The  absolute  [C  ii]  fluxes  are  among  the  faintest 
extragalactic  [C  ii]  emission  yet  seen.  In  Ligure  1,  we  plot  the 
[C  ii]  line  fluxes  against  the  LIR  flux  for  the  sample  ULIGs, 
Arp  299,  and  NGC  6240,  and  for  an  assortment  of  normal  and 
starburst  galaxies  presented  in  a  study  by  Lord  et  al.  (1996). 
The  LIR  flux  is  defined  as  described  above  for  both  the  ULIG 
and  normal/starburst  galaxy  samples. 

3.  DISCUSSION 

Observations  reveal  that  a  [C  n]-to-LIR  continuum  ratio  of 
0.1%-1%  characterizes  virtually  all  normal  and  starburst  gal¬ 
axies  (i.e.,  Crawford  et  al.  1985;  Stacey  et  al.  1991b),  as  il¬ 
lustrated  in  Ligure  1 .  Lor  normal  and  starburst  galaxies,  dense 


TABLE  2 

Sample  Properties  and  Measured  [C  ii]  Flux 


Source  Redshift  log  (LIR/L0)a  F„(60  (im)  F„(100  /xm)  Ffir  F[c, ]158„m 


Arp  299  .  0.104  11.91  103.7  107.4  4720  8.81  +  0.16 

IRAS  12071-0444  .  0.128  12.29  2.46  2.47  111  <0.05 

Mrk  231  .  0.042  12.54  31.99  30.29  1420  0.32  ±  0.04 

IRAS  15206+3342  .  0.124  12.18  1.77  1.89  81.2  0.20  ±  0.03 

IRAS  15250+3609  .  0.055  12.00  7.29  5.91  311  <0.14 

Arp  220  .  0.018  12.19  103.8  112.4  4790  0.87  +  0.04 

NGC  6240  .  0.024  11.82  22.68  27.78  1090  2.57  +  0.03 

IRAS  22491-1808  .  0.078  12.10  5.44  4.45  232  <0.12 

IRAS  23365+3604  .  0.065  12.13  7.09  8.36  336  0.122  ±  0.018 


Note. — F„(60  /xm)  and  Fv{  1 00  /xm)  are  from  IRAS  measurements  in  units  of  janskys.  FMR  and  F[c  II]15gjjm 
are  in  units  of  10~ 19  W  cm  "2.  Flux  upper  limits  are  3  a. 

“  Lir  defined  in  the  literature  as  the  IR  luminosity  from  8-1000  /xm;  data  from  Sanders  et  al.  1988a,  1988b, 
1991  and  Murphy  et  al.  1996. 
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Fig.  1. — The  [C  n]  158  fx m  line  flux  plotted  against  the  FIR  flux  for  the 
ULIGs  in  this  study  plus  Arp  299  and  NGC  6240  and  for  a  sample  of  normal/ 
starburst  galaxies  taken  from  Lord  et  al.  (1996).  The  error  bars  shown  denote 
the  [C  n]  flux  calibration  uncertainties,  which  are  much  greater  than  the  sta¬ 
tistical  errors.  The  two  dashed  lines  depict  the  observed  F[C  a]/FF1R  regime 
typical  of  normal  and  starburst  galaxies. 

(n  >  100  crrT3)  PDRs  exposed  to  far-UV  photons  from  OB 
stars  often  dominate  the  global  [C  ii]  luminosity  (e.g.,  Stacey 
et  al.  1991b).  Even  for  galaxies,  such  as  the  Milky  Way  and 
other  normal  spirals,  with  significant  [C  ii]  emission  from  the 
extended  low-density  warm  ionized  medium  and/or  the  cold 
neutral  medium,  F[Cii]/Ffir  still  ranges  from  0.1%  to  1%  (Mad¬ 
den  et  al.  1993;  Bennett  et  al.  1994;  Heiles  1994).  To  date, 
Malhotra  et  al.  (1997)  report  the  only  published  exceptions 
among  normal  and  starburst  galaxies  with  ratios  significantly 
below  0.1%,  where,  in  a  sample  of  30  normal  star-forming 
galaxies,  three  FIR-bright  sources  (LFIR~1010  L.-j  show 

L[Cii/Lfir<  2  x  10-4. 

Based  on  the  sources  in  this  study,  we  find  that  the  [C  ii] 
emission  relative  to  the  FIR  flux  in  UFIGs  is  generally  only 
10%  of  that  seen  in  normal  and  starburst  galaxies  (see  Fig.  1). 
For  six  of  the  seven  UFIGs  observed  to  date,  F[C  u]/Ffir  is  less 
than  5  x  10  4,  with  Arp  220  showing  the  lowest  measured 
ratio  among  our  sample  of  1.8  x  10  ’.  By  comparison,  Arp 
299  and  NGC  6240  display  [C  n]-to-FIR  continuum  ratios 
consistent  with  other  starburst  nuclei,  even  though  they  are 
only  slightly  less  luminous  mergers  than  the  UFIGs.  IRAS 
15206+3342  is  the  only  UFIG  presented  here  with  no  [C  n] 
deficit  (F[Cn]/FPIR  =  2.5  x  10  This  UFIG  is  a  Seyfert  2  gal¬ 
axy  exhibiting  “warm”  FIR  colors  of  F,,( 25  /xm )//-)/ 60  /im)  = 
0.2  and  F„( 60  /xm)/T'„(100  /tm)  =  0.9  (Sanders  et  al.  1988b); 
these  traits  are  not  unlike  those  of  other  UFIGs  in  our  sample, 
which  would  suggest  that,  within  the  UFIG  class,  the  [C  ii] 
deficit  does  not  depend  on  IR  luminosity,  FIR  color  temper¬ 
ature,  or  optical  spectral  type. 

In  the  ISO  FWS  FIR  spectra  of  Arp  220  and  Mrk  231,  the 
weak  [C  ii]  line  is  accompanied  by  other  weak  (undetected) 
fine-structure  lines  from  the  ionized  and  PDR  gas  (Fischer  et 
al.  1997a,  1997b,  1998).  These  lines  include  [O  m]  52,  88  /xm, 
[N  hi]  57  /im,  [N  ii]  122  /tm,  and  [O  i]  63,  146  /tm,  and,  like 
[C  ii],  are  generally  strong  in  the  FIR  spectra  of  starburst  gal¬ 
axies  such  as  Arp  299  (e.g.,  Satyapal  et  al.  1998).  A  comparison 
of  the  relative  strengths  of  lines  in  the  mid-  to  far-IR  spectra 
of  Arp  220  suggests  that  the  extinction  at  158  /tm  is  too  low 
to  explain  the  weak  [C  ii]  line  emission  (Fischer  et  al.  1997a, 
1997b,  1998).  With  this  in  mind,  we  consider  the  following 
possible  explanations  for  the  [C  n]  deficiency. 

Self-absorbed  or  optically  thick  [C  ii]. — Self-absorption  of 
[C  ii]  by  cooler  foreground  material  or  optically  thick  [C  n] 
emission  could  weaken  the  observed  [C  ii]  relative  to  FIR  in 


the  UFIGs.  In  Arp  220,  for  example,  Fischer  et  al.  (1997a, 
1997b)  report  self-absorption  of  the  [O  i]  63  /tm  line,  which 
is  a  strong  PDR  cooling  line  like  [C  ii].  In  our  Galaxy,  partial 
self-absorption  of  [C  ii]  has  been  recently  detected  toward  the 
molecular  cloud  associated  with  the  NGC  3576  H  ii  region 
(Boreiko  &  Betz  1997)  and  toward  SgrB2  using  the  ISO  FWS 
Fabry-Perot  interferometer  (Cox  et  al.  1998).  These  observa¬ 
tions  suggest  that,  particularly  for  edge-on  geometries,  [C  ii] 
self-absorption  could  be  important  in  UFIGs. 

Although  [C  ii]  is  generally  optically  thin  in  the  Galactic 
ISM  (e.g.,  Stacey  et  al.  1991a),  optically  thick  [C  ii],  i.e., 
tic  n]  158  >  2,  could  explain  the  observed  UFIG  F[C  u]/Ffir. 

Assuming  a  velocity  width  of  100-300  km  s  1  and 
NcJNh  =  1.4  x  10-4,  an  optical  depth  in  the  [C  ii]  line  of  2 
or  greater  is  reached  for  column  densities  NH  >  (2-6)  x  1023 
cnU2,  or  an  Av>  130-390  in  the  PDR  (C+)  gas  alone.  Such 
a  large  Av  disagrees  with  the  predictions  of  one-dimensional 
PDR  models  which  generally  show  that  the  [C  ii]  emission 
arises  from  a  thin  surface  layer  with  AVS4  (e.g.,  Tielens  & 
Hollenbach  1985).  Additionally,  recent  ISO  mid-IR  line  meas¬ 
urements  suggest  a  total  Av  of  ~50  toward  the  starburst  region 
in  Arp  220  (Sturm  et  al.  1996),  i.e.,  the  presumed  site  of  the 
[C  ii]  emission,  again  inconsistent  with  an  Av>  130-390  just 
in  the  PDR  zone. 

Saturation  of  [C  n]. — For  very  high  density  (n  »  /;cr  for 
[C  n]  =  3  x  103  cm  3)  PDRs  or  for  PDRs  with  a  high  ratio 
of  UV  flux  to  gas  density  (G0/n  +  1 0  cm3,  where  G0  is  nor¬ 
malized  to  the  solar  neighborhood  radiation  field;  Habing 
1968),  theoretical  PDR  models  show  that  7[ClI]  saturates  relative 
to  IFlR  such  that  /[C„//FiR  can  obtain  values  of  10~4  or  less  (e.g., 
Wolfire,  Tielens,  &  Hollenbach  1990;  Kaufman  et  al.  1998). 
In  the  high-density  case  where  n  lies  well  above  the  [C  ii] 
critical  density,  the  abundance  of  C+  in  the  upper  level  is  fixed 
by  the  atomic  statistical  weights.  In  addition,  the  column  den¬ 
sity  of  C+  becomes  insensitive  to  the  incident  UV  field  (Wolfire, 
Hollenbach,  &  Tielens  1989).  Since  both  the  abundance  and 
column  density  ceases  to  rise  with  n  and  G0,  the  emergent 
[C  ii]  line  intensity  saturates.  Within  a  PDR,  gas  heating  mainly 
occurs  via  the  ejection  of  energetic  electrons  from  dust  grains 
exposed  to  UV  photons.  In  the  high  GJn  regime,  the  high 
incident  UV  flux  increases  the  positive  charge  of  the  grain, 
thereby  reducing  the  energy  of  the  ejected  photoelectrons. 
Thus,  positively-charged  grains  reduce  the  fraction  of  UV  pho¬ 
ton  energy  that  is  converted  to  gas  heating,  which  causes  the 
[C  ii]  line  cooling  to  fall.  In  the  Galaxy,  [C  n]  line  saturation 
is  seen  toward  the  high-/;,  high-G()  PDRs  associated  with  Orion- 
like  H  ii  regions,  which  have  small  [C  n]-to-FIR  continuum 
ratios  of  about  a  few  times  1CT4  (Russell  et  al.  1980,  1981). 

Based  on  a  decreasing  trend  in  F[C  n]/FFm  with  increasing 
Ff  60  /im)//'’„(100  /xm),  i.e.,  increasing  rdust  or  G0,  Malhotra  et 
al.  (1997)  suggest  that  a  high  GJn  PDR  could  account  for  the 
[C  ii]  deficit  in  three  FIR-bright  galaxies.  For  the  UFIGs  ob¬ 
served  to  date,  we  see  no  such  trend  within  our  sample,  al¬ 
though,  as  with  the  Malhotra  et  al.  FIR-bright  galaxies,  the 
UFIGs  as  a  whole  are  characterized  by  low  F[CnJFmR  and  high 
(greater  than  0.8)  Ff  60  /im)//'’„(100  /im).  In  Arp  220,  we  com¬ 
pared  the  relative  strengths  of  the  [O  i]  146  /im  line  flux  upper 
limit  (Fischer  et  al.  1997b)  and  the  [C  ii]  and  FIR  fluxes  to 
the  predictions  of  the  PDR  models  of  Kaufman  et  al.  (1998) 
and  derived  two  solutions  for  n  and  G():  (1)  n  ~  200  cm-3, 
G0  ~  3  x  103  (GJn  ~  15  cm3)  and  (2)  n  ~  105  cm-3,  G()  ~  5 
(GJn  ~  5  x  10  5  cm3).  By  comparison,  other  normal  and  star- 
burst  galaxies  are  typically  characterized  by  G„/n  ~  0.1-1  cm3 
(e.g.,  Wolfire  et  al.  1990;  Carral  et  al.  1994;  Fischer  et  al.  1996), 
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which  suggests  that  if  [C  u]  line  saturation  alone  accounts  for 
the  low  F[C  n/FFIR,  then  the  PDR  physical  conditions  in  ULIGs 
differ  greatly  from  those  found  in  other  galaxies. 

Soft  UV  radiation  field. — A  stellar  population  deficient  in 
massive  main-sequence  stars,  either  due  to  an  aging  starburst 
or  an  initial  mass  function  with  a  low  upper  mass  cutoff,  could 
also  give  rise  to  a  soft  UV  field  and,  consequently,  a  small 
F[C  ii//V|r  (e.g.,  Satyapal  et  al.  1997).  In  this  case,  radiation 
from  cool  (Te{{  S  20,000  K)  stars  heats  the  dust  but  contains 
fewer  UV  photons.  Thus,  according  to  PDR  models  (Spaans 
et  al.  1994),  gas  heating  is  reduced  due  to  the  lower  rate  of 
photoelectrons  ejected  from  grains.  In  addition,  the  column 
density  of  C+  in  a  PDR  is  reduced  as  the  atomic/molecular 
transition  layer  moves  closer  to  the  cloud  surface.  The  net  effect 
reduces  the  size  and  luminosity  of  the  [C  ii]-emitting  region. 
If  the  [C  n]  deficit  results  from  an  age  effect,  then  there  must 
naturally  exist  ULIGs,  such  as  IRAS  15206+3342,  with  no 
[C  n]  deficit  that  presumably  contain  young  starbursts.  Like¬ 
wise,  an  unusually  high  dust-to-gas  ratio  could  diminish  the 
UV  flux  as  felt  by  the  gas,  since  the  dust  would  compete  more 
efficiently  with  the  gas  for  the  available  UV  photons.  This  effect 
would  decrease  the  extent  of  the  PDR  gas  and  thus  would  lower 
L[CiiJ/Lfir,  just  as  observations  of  low-metallicity  galaxies  show 
that  the  paucity  of  dust  can,  conversely,  increase  L[C  n]/LFIR  as 
a  result  of  the  greater  penetration  depth  of  the  UV  photons 
(e.g.,  Israel  et  al.  1996;  Madden  et  al.  1997).  Compared  to  a 
population  of  cool  stars,  a  high  abundance  of  dust  in  the  im¬ 
mediate  vicinity  of  the  starburst/ AGN  is  more  consistent  with 
the  warmer  LIR  colors  of  ULIGs,  although  some  unknown 
mechanism  for  increasing  the  dust-to-gas  ratio  is  required.  Un¬ 
like  the  aforementioned  explanations  for  the  [C  ii]  deficit,  a 
soft  or  severely-attenuated  UV  field  could  also  explain  the  lack 
of  LIR  [O  in]  lines  in  the  ISO  LWS  spectra  of  Arp  220  and 
Mrk  231  by  reducing  the  line  emission  from  the  ionized  gas 
as  well  (Lischer  et  al.  1997a,  1997b,  1998). 


4.  IMPLICATIONS 

As  examples  of  colliding  galaxies,  the  observed  ULIGs 
likely  represent  local  analogs  of  high-7  (£2)  protogalaxies.  If 
galaxies  form  from  the  coalescence  of  smaller  structures,  then 
protogalaxies  at  redshifts  from  1  to  10  should  resemble  the 
disturbed  ULIGs  at  low  z,  a  suggestion  that  is  supported  by 
the  appearance  of  the  faintest  sources  in  the  Hubble  Deep  Lield 
image  (e.g.,  Mobasher,  Rowan-Robinson,  &  Georgakakis  1996; 
van  den  Bergh  et  al.  1996).  As  the  brightest  line  in  the  spectrum 
of  most  galaxies,  the  [C  ii]  158  jitm  line  is  a  potentially  powerful 
spectroscopic  tracer  of  these  high-7  objects,  where  for  z  =  2-6, 
the  [C  ii]  line  is  shifted  to  submillimeter  wavelengths.  Recent 
derivations  of  the  [C  ii]  luminosity  in  protogalaxies  as  a  func¬ 
tion  of  z  suggest  that  it  may  be  possible  to  detect  as  many  as 
200  protogalaxies  per  winter  from  the  South  Pole  with  a  10  m 
telescope  (Stark  1997).  Such  an  estimate  assumes,  however, 
that  L[ClI]~  0.2%Lgal  from  observations  of  nearby  starburst  gal¬ 
axies.  The  much  smaller  L[CnJ/Lgal  ratio  implied  by  the  ULIG 
observations  presented  here  would  render  problematic  such 
ground-based  submillimeter  observations.  The  absence  of  sub¬ 
millimeter  [C  ii]  emission  from  the  z  =  4.69  QSO 
BR1033— 0327  would  seem  to  support  this  suggestion  (Isaak 
et  al.  1994).  Thus,  the  [C  n]-to-LIR  continuum  ratio  specific 
to  ULIGs  and  more  generally  to  protogalaxies  has  a  direct 
impact  on  the  feasibility  and  interpretation  of  future  submil¬ 
limeter  spectroscopic  surveys  from  the  ground  or  from  space. 
In  the  future,  we  will  combine  the  [C  n]  data  presented  in  this 
study  with  other  ISO  [C  ii]  observations  of  ULIGs  to  explore 
further  the  usefulness  of  [C  ii]  as  a  probe  of  high-7  galaxies 
and,  ultimately,  as  a  diagnostic  of  galaxy  formation  models. 
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